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We investigate the upper critical field in a stripe-phase and 
in the presence of a phenomenological pseudogap. Our results 
indicate that the formation of stripes affects the Landau orbits 
and results in an enhancement of H C 2- On the other hand, 
phenomenologically introduced pseudogap leads to a reduc- 
tion of the upper critical field. This effect is of particular 
importance when the magnitude of the gap is of the order of 
the superconducting transition temperature. We have found 
that a suppression of the upper critical field takes place also 
for the gap that originates from the charge-density waves. 
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I. INTRODUCTION 

The high-temperature superconductors (HTSC) ex- 
hibit qualitative differences with respect to the classical 
superconducting systems. Normal state properties of un- 
derdoped superconductors and the upward curvature of 
the upper critical field (H C 2 ) belong to one of the most 
spectacular examples. The presence of a normal-state 
pseudogap has been confirmed with the help of different 
experimental techniques like: angle-resolved photoemis- 
sion [jl] [|, intrinsic tunneling spectroscopy W5n, NMR 
H0, infrared || and transport || measurements. De- 
spite a wide spectrum of experimental data the underly- 
ing microscopic mechanism if far from being understood. 
A tempting hypothesis that the pseudogap is a precur- 
sor of the superconducting gap has not definitively been 
confirmed. In particular, the neutron scattering experi- 
ments ]l0| reveal qualitative differences between the iso- 
tope effects observed for the superconductivity and the 
pseudogap. Moreover, results obtained with the help of 
intrinsic tunneling spectroscopy ^,^] speak against the 
superconducting origin of the pseudogap. The coexis- 
tence of superconductivity and the charge-density waves 
pl|-|l4t can be considered as a possible scenario, that ac- 
counts for differences between superconductivity and the 
normal-state gap. Inhomogeneous distribution of holes 
which enter the copper-oxygen planes in the doping pro- 
cess can give rise to the formation of stripes. There is a 
convincing experimental and theoretical argumentation 
speaking in favor of the stripe-phase p5|-pl[| with an in- 
timate connection between superconducting and stripe 
correlations. This phase consists of antiferromagnetic 
domains which arc separated by the hole rich domain 
walls. 
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Differences between the high-temperature supercon- 
ductors and classical systems show also up in the mag- 
netic properties. The high-T c compounds are character- 
ized by large values of the upper critical field and its 
unusual temperature dependence. The resistivity mea- 
surements clearly indicate an upward curvature of the 
upper critical field with no evidence of saturation even 
at genuinely low temperatures |^,^3|. These results re- 
main in di sag reement with the conventional, microscopic 
approach |24| . This discrepancy can be explained as a re- 
sult of the Joscphson tunneling between superconducting 
clusters 1 25 2(| produced by a macroscopic phase separa- 
tion. 

Due to the complexity of the Gor'kov equations one 
usually assumes that the normal-state properties of the 
system under consideration can properly be described 
by three- |27| or two-dimensional 
cently, we have proposed an approac : 
culation of the upper critical field for a two-dimensional 
lattice gas ^9[[30) . This method allows one to derive 
in a similar way as one calculates the critical temperature 
in the standard BCS formalism. Therefore, any extension 
of the analysis of the upper critical field is rather straight- 
forward. In the present paper we calculate the upper 
critical field in a system that exhibits some important 
properties of hole-doped cuprates: stripe-phase and the 
presence of the pseudogap. In the latter case we discuss 
the coexistence of superconductivity and charge-density- 
wave as well as a phenomenological pseudogap. We show 
that the anomalous properties of the high-temperature 
superconductors are reflected in the upper critical field. 



II. H C 2 IN THE PRESENCE OF 
CHARGE-DENSITY-WAVES 

We consider a two-dimensional square lattice im- 
mersed in a perpendicular uniform magnetic field of mag- 
nitude H z . We assume the nearest-neighbor pairing 
interaction, Hy, that is responsible for anisotropic su- 
perconductivity and the interaction term, -ffcDW: which 
leads to the charge-density-waves. The relevant Hamil- 
tonian reads |E9,P0| 



where 



H — Hq + Hy + Hqdw , 
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Here, c^J annihilates (creates) an electron with spin a 
at the lattice site i, g is the gyromagnetic ratio and [1b 
is the Bohr magneton, (A) is the nearest-neighbor 
hopping integral that in the presence of the magnetic 
field acquires the Peierls phase-factor ]3l|j29| 



tij (A) = t exp 
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In the mean-field approximation the pairing interaction 
and the CDW coupling take on the form 
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where Ay = (cj^Cjf — c^Cji) is the superconducting sin- 
glet order parameter and <5cdw represents the magnitude 
of the CDW gap. The complexity of calculations strongly 
depends on the CDW modulation vector Q. For the sake 
of simplicity we consider a commensurate charge-density 
wave with Q = (jr, 7r). This choice of the modulation vec- 
tor results in the gap in the density of states that opens 
in the middle of the band (in our case at the Fermi level) 
independently on the magnitude of the external magnetic 
field. Since the pseudogap hardly depends on the mag- 
nitude of the mag netic field f§, <5cdw w iU be taken as a 
model parameter. 

In order to calculate the upper critical field we make 
use of the unitary transformation that diagonalizes the 
kinetic part of the Hamiltonian [ p9|j30| . In the case of 
the Landau gauge, A = H z (0, x, 0)7 this transformation 
is determined by a plane-wave function in y direction 
and an eigenfunction of the Harper equation p^] : 

g (k,p,m + l) + 2 cos (hm — pa) g (k,p, rrij 
+g(k,P,m- l) = t~ l E (k,p) g (k,p,m) . (6) 

Here, m is an integer number that enumerates the lattice 
sites in x direction, whereas h is the reduced magnetic 
field, h — 27r"I>/$0! that is expressed by the ratio of the 
flux $ through the lattice cell and the flux quantum <f>o. 
p is the wave-vector in y direction and k is an additional 
quantum number, that in the absence of the magnetic 
field is the wave-vector in x direction. In the new basis 
the normal-state Hamiltonian takes on the form: 



^cdw = ~(5cdw E g* (k,p + ir,m) g (l,p, 
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where 
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One can prove that if E represents the eigenvalue of the 
Harper equation obtained for the wave-vector p then — E 
is one of the eigenvalues corresponding to p + n. It can 
also be shown that g (k, p, rnj = g (k, p + tt, m) exp (?7rm) 
represents an eigenfunction of the Harper equation calcu- 
lated for momentum p with the eigenvalue — E (k, p + 7r) . 
With the help of these relations one can obtain analyti- 
cally the energy spectrum of the Hamiltonian in the nor- 
mal state, Hq + -ffcDWj provided that the eigenvalues 
of the Harper equation are known. In particular, one 
can calculate the anomalous Green functions which are 
related to the superconducting order parameter: 
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where X and Y are determined by the solution of the 
Harper equation 

x Ik P ( m ) = 9 (I m) g (k, -p, m + l) 

+g(l,p,m + l) g(k,-p,m) , (11) 
Y lk P ( m ) = 2cos(p)g (l,p, m) g (k, -p, m) . (12) 

Due to the plane-wave behavior in y direction the super- 
conducting order parameter depends explicitly only on 
the position in x direction (see Ref. |Q for the details) 
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Kjj. p (u)), when integrated over w with the Fermi function, 
gives the Cooper pair susceptibility 



u + gjgt) fa - E k,-pl) + ^CDW 
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(15) 

One can see that the impact of the charge-density waves 
on superconductivity is brought about only by the mod- 
ification of this quantity. Equations (10), (13) and (14) 
allow one to calculate the upper critical field. It is de- 
termined as the highest magnitude of the magnetic field 
for which there exists a non-zero solution for A x (m) and 
Av(m). 
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FIG. 1. Temperature dependence of the upper critical field 
for different magnitudes of the CDW gap. 

Figure 1. shows the reduced upper critical field, h c2 , 
as a function of temperature calculated for different mag- 
nitudes of the CDW order parameter, <5cdw- We have 
adjusted the magnitude of the pairing potential V that 
gives the same superconducting transition temperature 
kT c = 0.02i for all values of <5cdw- These results have 
been obtained for 120 x 120 cluster that at temperatures 
kT ~ 10~ 2 i gives convergent results (we refer to Ref. [3(J 
for details of the cluster calculations). One can see that 
even for small magnitudes of the CDW order parameter 
the upper critical field is significantly reduced. However, 
qualitative temperature dependence of H C 2 is not affected 
by the charge-density wave correlations. The reduction 
of the upper critical field due to the charge-density waves 
may be brought about by a direct coupling between CDW 
and superconducting order parameters as well as by the 
modification of the density of states. In order to distin- 
guish these contributions we investigate the upper critical 
field in the presence of a phenomenological normal-state 
gap of arbitrary magnitude and depth. This problem will 
be discussed in the next section. 



III. H C 2 IN THE PRESENCE OF A 
PHENOMENOLOGICAL GAP 

In this section we investigate modification of the upper 
critical field that originates only from the normal-state 
gap in the density of states. In contradistinction to the 
analysis presented in the previous section, the density 
of states may remain finite despite the presence of the 
gap. Here, the normal-state gap is characterized by the 
width 26 and the relative depth (po — Ppg) I Pa> where 
Ppg and po denote the density of states in the presence 
and without the pseudogap, respectively. It is visualized 
in the inset in Fig. 2. In the absence of the CDW order 



the upper critical field is determined by Eq. (10) with 
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In order to account for the modification of the density 
of states we renormalize the normal-state propagators 
which give rise to the Cooper-pair susceptibility 
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In the limiting case ppc = Po one obtains the standard 
density of states as determined by the Hofstadter spec- 
trum, whereas for ppc = the density of states van- 
ishes in the vicinity of the Fermi energy. Substituting the 
renormalized propagators into Eq. (16) one can calculate 
the upper critical field in the same way as described in 
the previous section. 
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FIG. 2. Critical temperature as a function of the pseudogap 
width S for different values of the external magnetic field. 
Ppg = §Po was used. The inset shows a schematic density of 
states in the vicinity of the Fermi level. 

Figure 2. shows the superconducting transition tem- 
perature obtained for different values of the reduced mag- 
netic field with ppg = \po- As before, the intersite cou- 
pling V have been adjusted to obtain kT c = 0.02t in 
the absence of magnetic field. One can see that the up- 
per critical field is reduced due to the presence of the 
normal-state gap. The most significant lowering of H c2 
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takes place for finite values of the S which are compara- 
ble to the magnitude of the superconducting gap. This 
result originates from the fact that the Cooper pair sus- 
ceptibility is strongly peaked at the Fermi level with a 
characteristic energy scale that is determined by temper- 
ature. Therefore, for S » kT c the pseudogap results in 
a global lowering of the density of states which can be 
compensated by an enhancement of the pairing poten- 
tial. It means that assuming stronger pairing potential 
V we can reproduce H C 2(T) calculated in the absence of 
the gap. 
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FIG. 3. Upper critical field as a function of temperature 
for different densities of states at the Fermi level. The half 
width of the pseudogap is 5 — O.lt. The inset shows tem- 
perature dependence of upper critical field calculated for the 
phenomenological pseudogap with ppG = (solid line) and 
for the CDW gap (dashed line). In both cases the half width 
is S — Scow = O.Olt. 

We have found that the reduction of the upper critical 
field increases with the depth of the gap as depicted in 
Fig. 3. The inset in Fig. 3 shows a comparison of H C 2 ob- 
tained for the phenomenological pseudogap with ppq = C 
and for the charge-density waves. A comparison of these 
results clearly indicate that the coupling between CDW 
and superconducting order parameters results in a small 
decrement of the upper critical field. 



IV. H C 2 IN A STRIPE PHASE 

Other unusual feature of HTSC, that we discuss in 
the present section, is related to inhomogeneous distri- 
bution of holes. It results in a stripe-phase which con- 
sists of antifcrromagnetic domains separated by hole-rich 
domain walls. We study how the upper critical field is 
affected by this specific distribution of carriers. In order 
to simulate the presence of a stripe-phase we carry out 



the calculations for a long and narrow rectangular-shape 
clusters. We assume that the isolating, antifcrromagnetic 
domains can be simulated by fixed boundary conditions 
in the direction perpendicular to the stripes (along the 
x axis). The spatial organization of the stripe structure 
has been intensively investigated on experimental [j33|34 
and theoretical grounds |35 36 1. Experimental data for 
HTSC show that the width of stripes depends on the 
concentration of holes and is of the order of a few lattice 
constants. The neutron-scattering study of the stripe 
phase (3^] suggests that the hole-rich domain walls are 
only single cell wide. On the other hand, the numeri- 
cal study of the two-dimensional t-J model |55| shows 
that the domain walls may have a significant density of 
holes over three rows of sites. According to these results 
we consider 150 x n finite systems, where n = 2,3 and 
7. Our simplified approach does not restore the actual 
structure of the stripe-phase. In particular, for n = 1 
one obtains an unphysical, purely one-dimensional sys- 
tem, that hardly depends on the external magnetic field. 
Therefore, we investigate the rectangular-shape clusters 
with the width as a free parameter. Since we neglect the 
correlations between different stripes, the upper critical 
field is determined by Eqs. (10) and (16). 

In the case of free electron gas external magnetic field 
leads to the occurrence of rotationally invariant states 
corresponding to the Landau orbits. However, the ge- 
ometry of the stripe-phase may seriously affect the for- 
mation of the Landau orbits. This effect is of particular 
significance if the radii of the Landau orbits, Rl, exceed 
the width of the stripe, an, (a is the lattice constant). In 
order to visualize the impact of magnetic field on elec- 
trons in the stripe-phase we have calculated the result- 
ing current distribution. Within the framework of the 
linear-response theory the current operator is given by 
Ji(x,y) = — dH / dAi{x, y), where {x,y) denotes spatial 
coordinates and Z's are unit vectors in the lattice axes 
directions. Results obtained in the normal state (V = 0) 
on a 150 x 7 cluster with applied magnetic field h = 0.1 
are presented in Figure 4. 
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FIG. 4. Current distribution in a piece of stripe of width 
of 7 sites, calculated for the reduced magnetic field h = 0.1. 
The lengths of the arrows are proportional to the currents. 
Such a pattern is periodically repeated along the stripe. 
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Modification of the Landau orbits affects the diamag- 
netic pair-breaking mechanism. Therefore, one may ex- 
pect that superconductivity survives in the presence of 
much stronger magnetic fields than in the homogeneous 
phase. This observation is confirmed by the numerical 
calculations, as depicted in Figure 5. 




FIG. 5. Upper critical field as a function of temperature 
calculated for stripes of different width. We have chosen ap- 
propriate values of V which give the same transition tem- 
perature in the absence of magnetic field. The inset shows a 
comparison of upper critical field for 150 x 150 and 150 x 7 sys- 
tems with the same value of the pairing strength V = 0.244t 

Here, 150 x 150 cluster corresponds to an infinite system. 
The enhancement of H C 2 is of particular importance for 
weak magnetic fields, when R^/na — > oo. One can ob- 
serve a dramatic change of the slope, dTJ^/dT, calcu- 
lated at T — T c . Here, the impact of the magnetic field 
on the superconducting transition temperature is much 
less than in the homogeneous two-dimensional case. 

The pseudogap and stripes affect the superconduct- 
ing properties of the system both in the presence and in 
the absence of the magnetic field. Modification of the 
density of states changes the effective coupling constant, 
A = Pfs^, that enters the standard BCS gap equation. 
Therefore, we have directly compared the H C 2 for sys- 
tems, which in the absence of magnetic field are char- 
acterized by the same transition temperature (one can 
roughly say that A = const). In order to complete the 
discussion, we have also calculated the H C 2(T) for the 
case when the pairing potential does not depend on the 
pseudogap and the stripe structure (V = const). Since, 
the opening of the pseudogap reduces T c it results also in 
an additional decrement of the upper critical field, when 
compared to the results presented in Figs. (1-3). How- 
ever, an enhancement of the H C 2 in the stripe phase can 
take place despite the reduction of the superconducting 
transition temperature, as depicted in the inset in Fig. 
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V. CONCLUDING REMARKS 

In order to clarify some physical aspects of our method 
one can compare it with approaches, which are commonly 
used to investigate H c %. Previously, we have applied 
the same method to discuss the upper critical field for 
isotropic superconductivity p9| , p0| . Then, one ends up 
with the gap equation that can be written in the form 

A, = j AjG(i, j>„)G(z, j, -cu n ). (18) 

Here, A^ = (ciic^) and G(i,j,uj n ) is the one-electron 
Green's function in the presence of a uniform and static 
magnetic field. It is clear that the above equation is a 
lattice version of the linearized Gor'kov equations 
which determine the critical field at a second-order tran- 
sition, where the superconducting gap A^ is vanishing 
. However, our method does not allow to discuss the 
superconducting properties below the H C 2 (e.g. the vor- 
tex state). In our approach the electron Green's functions 
have been calculated exactly, whereas in the standard 
case one makes use of the semiclassical approximation 
that neglects the Landau level quantization. 

To conclude, we have investigated how the upper crit- 
ical field is connected with different features of high- 
temperature superconductors. In particular, we have 
discussed H c2 in the presence of charge-density waves, 
phenomenological pseudogap and stripes. Our results 
suggest that a gap in the density of states reduces the 
upper critical field, independently on the underlying mi- 
croscopic mechanism. For finite density of states at the 
Fermi level this reduction is mostly pronounced when 
the width of the gap is of the order of the superconduct- 
ing transition temperature. In the phase with isotropic 
CDW gap the density of states at the Fermi level van- 
ishes. Then, as one can expect, the upper critical field is 
strongly reduced even by a relatively small gap. Here, the 
coupling between the CDW and superconducting order 
parameters results in an additional reduction of H C 2- On 
the other hand, in the presence of stripes the upper criti- 
cal field is enhanced, especially close to T c . We attribute 
this effect to the reduction of the orbital pair-breaking 
mechanism since the radii of the Landau orbits are much 
larger than the width of the stripes. 

The presented investigation of H C 2 is restricted to the 
simplest case of the uniform magnetic field and neglects 
a possible disorder in the vortex system. It can originate 
from fluctuations close to the phase transition or inhomo- 
geneous charge and spin distribution in the stripe-phase. 
However, as we are concerned exclusively with the criti- 
cal field at the second-order transition, these effects are 
of minor importance. We have also not discussed the 
reentrance of superconductivity in the strong-magnetic 
field. This effect has been investigated in the continuum 
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model 38 3^] as well as in the case of lattice gas [ p9| , 
when the structure of fractal energy spectrum is reflected 
in phase diagram. Theoretical argumentation that sup- 
ports the reentrance of superconductivity remains valid 
also in the presence of pseudogap, at least on the simplest 
level that has been used in the present paper. However, 
in the genuinely strong magnetic field the assumption of 
the field-independent gap is unphysical and microscopic 
investigation of this phenomenon is needed. 
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